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^H Abstract 

O 

CN To date, Bi2Se3 is known as the best three-dimensional high temperature topological insulator with a 

^ large bulk band gap. Here, we theoretically predict that the series of Pb-based layered chalcogenides, 

O Pb^Bi2Se^+3 and Pb^Sb2Te^+3, are possible new candidates for topological insulators. As n increases, 

_ the phase transition from a topological insulator to a band insulator is found to occur between n = 2 

and 3 for both series. Significantly, among the new topological insulators, we found a bulk band gap of 

■J^ 0.40eV in PbBi2Se4 which is one of the largest gap topological insulators, and that Pb2Sb2Te5 is located 

C in the immediate vicinity of the topological phase boundary, making its topological phase easily tunable 

53 by changing external parameters such as lattice constants. Due to the three-dimensional Dirac cone at the 

S 

phase boundary, massless Dirac fermions also may be easily accessible in Pb2Sb2Te5. 
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Topological insulators, distinguished from normal band insulators by a nontrivial Z2 topolog- 
ical number and topologically protected surface states, have attracted great attention due to their 
significance both for applications and for fundamental research on a new quantum state of matter 
|[Tl[2l]. Since the suggestion of a quantum spin Hall effect in a honeycomb lattice EIH, which is 
a time-reversal pair of Haldane models iH and induced by spin-orbit interactions, several topo- 
logical insulator materials were predicted and observed [l6l-[T2ll. Efforts to find new topological 
insulators are still being made [[I3l[l4l. A large spin-orbit coupling (SOC) strength is essential to 
realizing non- trivial topological band structures. 

To classify the Z2 topological number, we adopted the method proposed in ref. ifTSll . which 
studied the topological phase transition by making a variation of external parameters such as the 
SOC strength or lattice constants. In our work, we varied the SOC strength (Aso) from to the real 
values of the systems (Aq), and investigated whether the topological phase transition occurs. In 
three-dimensional (3D) topological insulator materials, a band crossing between conduction and 
valence bands should occur during this process: the 3D Dirac cone appears at the transition point 
between the band insulator (BI) and topological insulator (TI). 

There are several advantages of this method. First, it is intuitive and only needs a direct obser- 
vation of the Dirac cone at the transition point. Second, it does not require a heavy computational 
cost because of the simple unit-cell calculations. The surface calculations to see the topologically 
protected surface state are also one of alternative ways to determine the topological phase. The 
degrees of complexity of such calculations are, however, grater than the method we employed 
here. Third, this method can be applied to every system regardless of the existence of inversion 
symmetry. Lastly, the critical SOC value at which the band crossing occurs shows how far the 
system is located from the phase boundary of BI and TI. 

By using this method, we found that the series of Pb-based chalcogenides, Pb^Bi2Se^+3 and 
Pb^Sb2Te^+3, are possible new candidates for TI materials and the topological phases are changed 
from TI to BI with increasing n. Among the above series, we focus on the material near the phase 
boundary between BI and TI due to the possibility of tuning the topological phase by changing 
external parameters. 

To investigate the electronic structures and topological phases, first-principles calculations were 
performed using the full-potential linearized augmented plane wave method [[T6l with the gradient- 
corrected Perdew, Burke, and Ernzerhof form of the exchange-correlation functional [17J. The 
core states and the valence states were treated fully relativistically and scalar relativistically, re- 



spectively. The calculations were carried out with the experimental lattice parameters for bulk 
PbBi2Se4, PbSb2Te4 and Pb2Bi2Se5, and with the fully optimized geometry for the others, because 
crystal structures of n = 1 compounds and Pb2Bi2Se5 were reported experimentally [fT8] - [20ll and 
those of other n's are designed here theoretically. 

PbBi2Se4 and PbSb2Te4 have the rhombohedral crystal structure, and a layered structure 
stacked along the c-axis of the hexagonal lattice, consisting of seven atoms in one septuple layer. 
The Pb atom is sandwiched by Se-Bi-Se or Te-Sb-Te layers and located at the inversion center, 
shown in Fig. 1(a). There are van der Waals interactions between two septuple layers, and this 
provides a natural surface geometry which is appropriate for observing a topologically protected 
surface state. 

Upon increasing the SOC strength ratio, Aso/^o, from to 1, both materials show the phase 
transition from BI to TI; in other words, there are band crossings during the process, cf., Figl.(b)- 
(g). Critical values of the ratio are 0.46 for PbBi2Se4 and 0.62 for PbSb2Te4, and the electronic 
band structures at each critical point are shown in Fig. 1(c) and Fig. 1(f). 3D Dirac cones are seen 
at the Z point, one of the time-reversal invariant momentum points where -k is equivalent to k; 
they are doubly degenerate due to spatial inversion and time reversal symmetry. 

Without the spin-orbit interaction, the system should have a trivial Z2 topological number (u = 
0), the so called BI. The presence of band crossings between conduction and valence bands during 
increasing SOC represents the change of its Z2 topological number from trivial (u = 0) to non- 
trivial (u = 1) insulators. If the topological phase transition occurs before the SOC strength 
reaches Aq, the system is classified as TI. In these inversion symmetric systems, band crossing is 
equivalent to band inversion where conduction and valence bands exchange their parities at the 
crossing point. 

From our calculations, the bulk gaps of PbBi2Se4 and PbSb2Te4 are 0.28eV and 0.13eV, respec- 
tively. The gap size of PbBi2Se4 is comparable to that of Bi2Se3 [El [IS. The gap values from 
the calculations without SOC are 0.24eV (Fig. 1(b)) and 0.26eV (Fig. 1(e)), which reflects purely 
hybridization effects. After increasing the SOC strength, gaps shrink to zero and then increase 
again, which shows the competition between hybridization and SOC strength in the topological 
phase transition. A non-trivial topological phase is a result of the predominance of SOC over 
the hybridization strength. Compared to other TI chalcogenides, the large band gap in PbBi2Se4 
originates from weak hybridization and large SOC strength. 

To estimate the gap size of PbBi2Se4 more precisely, we performed non-local screened- 



exchange LDA (sX-LDA) calculations [[2T1 |22l which overcome the well-known weakness of 
gap underestimation in the LDA scheme and reproduces the experimental gap values in many 
semiconducting materials. Significantly, the sX-LDA result yields a gap of 0.40eV in PbBi2Se4, 
which confirms that PbBi2Se4 is one of the largest gap TI, leading to the fact that PbBi2Se4 may 
be the best for high temperature applications among the known TI materials. 

One advantage of the method, used in this work to verify the topological phase, is that it de- 
scribes the exact distance from the phase boundary of BI and TI in terms of the SOC strength. For 
instance, the larger critical SOC strength of PbSb2Te4 than that of PbBi2Se4 means that the former 
is closer to the phase boundary than the latter. 

In addition to PbBi2Se4 and PbSb2Te4, we suggest the series of Pb^Bi2Se^+3 and Pb^Sb2Te^+3 
structures where n is an integer larger than 1 and show that a phase transition from TI to BI 
occurs as n increases. The crystal structure of Pb^Bi2Se^+3/ Pb^Sb2Te^+3 is composed of 
(PbSe)^/(PbTe)^ core with Se-Bi-(Se)/Te-Sb-(Te) sandwich layers. As shown in Fig.2(a), the 
critical SOC ratio (Ago/Ao) exceeds 1 f or n > 3 for both series of materials, which means that 
there is a phase boundary between n = 2 (TI) and n = 3 (BI). The gap sizes decrease up to n = 2 
and increase again in going beyond the phase boundary, (cf. Fig. 2(b)) On the other hand, the 
calculated gap values without SOC rise monotonically with n; this shows that as n increases, the 
hybridization strength exceeds SOC, resulting in the topological phase transition. 

In Fig. 2(c) and (d), electronic band structures from the slab geometry of the n = 1 and n = 3 
composition of Pb^Bi2Se^+3, which consist of 6 septuple and 4 hendecuple layers with 78.4A and 
85. SA respectively, are shown. Consistent with the critical SOC ratio in Fig. 2(a), PbBi2Se4 has 
a topologically protected surface state, whereas there is no such state connecting the valence and 
conduction bands in the Pb3Bi2Se6 slab. This is further evidence of the change in Z2 topological 
number from 1 to as n increases. The 2D surface state Dirac cone in the PbBi2Se4 slab is robust 
under the presence of non-magnetic perturbations. Again, the 0.35eV gap originating from bulk 
states is useful for high temperature spintronics applications. In other words, a single isolated 
surface state Dirac cone is seen in the energy range from -O.OSeV below the Fermi level to 0.32eV 
above. In the case of the Pb3Bi2Se6 slab, due to the weak van der Waals type interlayer interactions, 
no distinct surface state appears from the bulk BI states. 

Among the Pb-based chalcogenides series, Pb2Sb2Te5 whose critical SOC ratio is 0.96, is lo- 
cated quite close to the phase boundary, that its Z2 topological number may be easily tuned by 
changing external parameters. The 26meV gap in the fully optimized geometry of Pb2Sb2Te5 



(Fig. 3(a)) is closed by making a small change in lattice constants. As shown in Fig. 3(b), Pb2Sb2Te5 
with a 0.5% reduction of the a6-lattice constants is positioned exactly at the borderline of TI 
and BI, whose electronic band structure shows a gapless 3D Dirac cone at the F point, realizing 
monopoles which generate an inverse- square type Berry curvature in the 3D Brillioun zone. A re- 
duction of the a6-lattice constants makes the hybridization stronger, resulting in equilibration with 
the SOC strength. The band structure near the Dirac point shows anisotropic dispersions between 
in-plane and out-of-plane directions, described by a massless Dirac Hamiltonian 

H = v±a± • k± + '^||(j|| • fcy (1) 

where a are Pauli matrices, k± is a momentum vector along the c-axis and k\\ is on the a6-plane. 
Nearly circular cross-section on the (kx-ky) plane and an ellipsoidal cross-sections on the (k^-kx) 
plane are depicted in Fig. 3(c) and (d), where the electronic band dispersion and the iso-energy 
contours of the 3D Dirac cone are shown. 

The substitution of Te by Se might be useful to control the topological phase of this material 
in terms of varying both the SOC strength and lattice parameters. Considering the smaller ionic 
radius and SOC strength of Se than those of Te, the Se substitution can play a role in reducing 
not only lattice constants but also the SOC strength which can both accelerate the phase transition 
from TI to BI. Hence, Pb2Sb2(Tei_a;Sea;)5 might be appropriate for studying the topological phase 
transition and the 3D Dirac cone. 

In this work, we predicted new TI materials of Pb-based chalcogenides by investigating the 
change of the Z2 topological number during variation of the SOC strength. With increasing n 
in the Pb^Bi2Se^+3 and Pb^Sb2Te^+3 series, the TI to BI transition occurs between n = 2 and 
3. Among these compounds, Pb2Sb2Te5 has a critical ratio of SOC strength close to 1, which 
means closeness to the phase boundary. Both topological phases and the gap at the 3D Dirac 
cone are tunable by changing lattice constants or the SOC strength. Also, it will be interesting 
to investigate the physical properties of the bulk 3D Dirac cone manifested at the critical point. 
As an extension of the present work, the Pb^(Bi/Sb)2m(Se/Te)^+3^ series, which are pseudo- 
binary [Pb(Se/Te)]^[(Bi/Sb)2(Se/Te)3]^ systems and combinations of Pb^(Bi/Sb)2(Se/Te)^+3 and 
[(Bi/Sb)2(Se/Te)3]^ layered structures, are potential candidates for new TI materials which might 
show various bulk band gaps and topological phase transitions. The method used in our work 
can be the most appropriate tool to determine the Z2 topological number of the systems without 
inversion symmetry, which show pair creation and annihilation of monopole and anti-monopoles 



at the topological phase transition region lITSll . 
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FIG. 1. (a) Crystal structures of PbBi2Se4/PbSb2Te4. Electronic band structures of (b)-(d) PbBi2Se4 and 
(e)-(f) PbSb2Te4. Calculations (b),(e) without spin-orbit interactions (Aso = 0), (d),(g) with real SOC 
strength (Aso = Aq), and (c),(f) with critical SOC strength (Aso = -^so)- 
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FIG. 2. (a) Critical SOC ratio and (b) bulk band gap with respect to n for the Pb„Bi2Se„+3 and 
Pb„Sb2Ten+3 series. Band structures of slab geometry of (c) 6 septuple layers of PbBi2Se4 and (d) 4 
hendecuple layers of Pb3Bi2Se6, which represent TI and BI, respectively. 
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FIG. 3. Band structures of Pb2Sb2Te5 with (a) fully optimized crystal structure and (b) 0.5% reduction of 
a6-lattice constants. Three-dimensional plot of Dirac cone in (b) and its contours are shown in (c) and (d). 
Anisotropic Dirac cone dispersions on (c) kx-ky and (d) kz-kx plane. 



